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Oh. Wc prove H^ orbital stability of Dirac solitons in the integrable massive Thirring model by 

■^^ ' working with an additional conserved quantity which complements Hamiltonian, momentum 

ly-v . and charge functionals of the general nonlinear Dirac equations. We also derive a global bound 

on the H^ norm of the L^-small solutions of the massive Thirring model. 

<^ ; 1 Introduction 



Nonlinear Dirac equations are considered as a relativistic version of the nonlinear Schrodinger 

(NLS) equation. Coinpared to the NLS equation, proofs of global existence and orbital stability 

of solitary waves are complicated by the fact that the quadratic part of the Hamiltonian of the 

nonlinear Dirac equations is not bounded from neither above nor below. Similar situation occurs 

^1 in a gap between two bands of continuous spectrum in the Schrodinger equations with a periodic 

OO ' potential, for which the nonlinear Dirac equations are justified rigorously as an asymptotic model 

^ : (Chapter 2.2 in [35]). 

Because orbital stability of solitary waves is not achieved by the standard energy arguments jl7] . 
researchers have studied spectral and asymptotic stability of solitary waves in many details. Spectral 
f^ . properties of linearized Dirac operators were studied by a combination of analytical methods and 

^^ \ numerical approximations |15| [H [3l [6l [im [TTl 112] . Asymptotic stability of small solitary waves in 

the general nonlinear Dirac equations was studied with dispersive estimates both in the space of 
one [Ml [261 123 and three [H [3 [7] dimensions. Global existence and scattering to zero for small 
^ \ initial data were obtained again in one [TH [19] and three [29l [30] dimensions. 

When nonlinear Dirac equations are considered in one spatial dimension, a particular attention 
is drawn to the massive Thirring model (MTM) [38], which is known to be integrable with the 
inverse scattering transform method [25^ I28j . In laboratory coordinates, this model takes the 
following form: 

i{ut + Ux) + v = 2\v\'^u, ,. 

i{vt — Vx) -\-u = 2\u\'^v, 

where {u,v){x,t) : M x M+ ^ C^. 

Selberg and Tesfahun [37] proved local well-posedness of the MTM system in H^iM.) for s > 
and global well-posedness in i?*(R) for s > 2- Machihara et al. [31] proved for similar nonlinear 
Dirac equations with quadratic nonlinear terms that local well-posedness holds in H'^ (M) for s > — 2 
and that the Cauchy problem is ill-posed in i?~^'^(]R). Candy [8] proved local and global well- 
posedness of the MTM system in L^(M). These results do not rely on the inverse scattering 
transform for the MTM system. 



On the other hand, by using the inverse scattering transform, spectral stabihty of the MTM 
sohtons was estabhshed by Kaup & Lakoba [23]. The stationary MTM sohtons are known in the 
exact analytical form: 



^; = [/^(x + 2;o)e*^*+*°, ^^^ 



with 



vr 



(jj 



2 



Uu^ix) = / N ; / N ^ (3) 

VT+tJcosh ( V 1 — uj'^x J + i^/l^^uJ sinh. f y 1 — uj'^x I 

where a and xq are real parameters related to the gauge and space translations, whereas cj G (—1,1) 
is a parameter that determines where the MTM sohtons are placed in the gap between two branches 
of the continuous spectrum of the linearized MTM system ([1]) . Perturbations of the MTM system 
and the loss of spectral stability of solitary waves was consequently considered using the spectral 
representations [21 [2l] and the Evans function [21] , No results on the orbital or asymptotic stability 
of the MTM solitons ([2]) in the time evolution of the MTM system ([T]) have been obtained so far. 
The idea for our work relies on the existence of an infinite set of conserved quantities in the MTM 
system, which has been known for quite some time [28]. The three standard conserved quantities 
for the nonlinear Dirac equations are related to the translational invariance of the system with 
respect to gauge, space, and time transformations. For the MTM system ((Tj), these three conserved 
quantities are referred to as the charge Q, momentum P, and Hamiltonian H functionals: 

Q= f (|n|2 + |^;|2) dx, (4) 



and 



P = - (uux - UxU + vvx - Vxv) dx, (5) 

^ Jm. 

H = - I {uux — UxU — vvx + Vxv) dx + [—vu — uv + 2\u\ \v\ ) dx. (6) 

2 Jr Jr 

The charge Q is useful to establish global bound on the L^ norm of solutions, as soon as the local 
existence in L'^{M.) is proven [8]. The other two functionals P and H are defined in H^''^{M.) but 
they are not so useful because of the fact that the quadratic part of the Hamiltonian H is not 
bounded from neither above nor below. Therefore, although the MTM solitons ([2]) are critical 
points of the functional H + uQ + cP, where w E (— 1, 1) is a parameter of the MTM solitons ([3|) 
and c = for stationary solutions, this functional cannot serve as a Lyapunov functional for orbital 
stability or instability of the MTM solitons. 

Nevertheless, we find another conserved quantity of the MTM system ([T|), thanks to the inte- 
grability via the inverse scattering transform method: 



R 



\ux\'^ + |wxP - -{uxU - Uxu){\u\'^ + 2|u|2) + -{vxV - Vxv){2\u\'^ + l^p) 
-{uv + uv){\u\^ + |7;p) + 2|'up|z;p(|n|2 + \v\^)] dx. (7) 



Derivation of the conserved quantity R is reviewed in Appendix A. The conserved quantity R is 
well defined in H^(R) and we shall use it to prove orbital stability of the MTM solitons in H^(R). 
The main result of this article is the following theorem. 



Theorem 1. There is coq G (0, 1] such that for any fixed uj G (— wq, wq), the MTM soliton {u,v) = 
{Uuj, U^) is a local non- degenerate minimizer of R in i7^(]R, C^) under the constraints of fixed values 
of Q and P. Therefore, the MTM soliton is orbitally stable in H^(R,C'^) with respect to the time 
evolution of the MTM system U^). 

Prom a technical point, we estabUsh that the MTM sohtons ([2]) are critical points of the func- 
tional Au) = R+ {1 — uj'^)Q, where to G (—1,1) is the same parameter of the MTM solitons ([3]). By 
using operator calculus in constrained spaces, we prove that there is ujq G (0, 1] such that for any 
fixed CO G {—coq, loq), the functional A^j is strictly convex at the MTM soliton (n, v) = {U^^, Uuj) under 
the constraints of fixed values of Q and P. As a result, A^^ can serve as a Lyapunov functional for 
orbital stability of the MTM solitons thanks to the conservation of R, Q, and P and the standard 
analysis of orbital stability of solitary waves [17]. Appendix B states relevant results used in our 
work. 

Note that the non-degenerate minimizer in Theorem [1] can be translated along two "trivial" 
parameters a and xq in ([2]). These parameters are related to the gauge and space translations 
and can be excluded by additional constraints on the perturbations to the MTM soliton (u, v) = 
{Uuj, tJui)- Whereas we have not succeeded to find the exact value of wq, we conjecture that ujq = 1, 
that is, the result of Theorem [T] extends to the entire family of MTM solitons. 

We also mention some recent relevant results. First, orbital stability of breathers of the modified 
KdV equation is proved in space H'^{M) in the recent work of Alejo and Munoz ^. An additional 
conserved quantity is introduced and used to complement conservation of the Hamiltonian and 
momentum of the modified KdV equation. This work is conceptually similar to the ideas of our 
paper with the following difference. It uses the known characterization of orbital stability of multi- 
solitons in the KdV and modified KdV equations with higher-order conserved quantities \1Q\ [32] , 
whereas our work introduces a new concept of orbital stability of Dirac solitons. 

Second, a different technique involving additional conserved quantities is proposed in the work of 
Deconinck and Kapitula [H], where no constraints are imposed to study orbital stability of periodic 
waves with respect to perturbations of multiple period in the KdV equation. The lack of minimizing 
properties of the higher-order Hamiltonian is corrected by adding lower-order Hamiltonians with 
specially selected strength parameter. 

As a bi-product of our work, we obtain a global apriori bound on L^-small solutions of the 
MTM system in H^{E.). The standard apriori bounds on the H^ norm of the solution of a general 
nonlinear Dirac equation grow at a double-exponential rate [161 134] . At the present time, we do 
not know if global bounds on the H^ norm can be proven for all (not L^-small) solutions of the 
MTM system ^. 

We also do not know if scattering to the MTM solitons (and hence asymptotic stability of the 
MTM solitons) can be proved by using the inverse scattering transform methods, e.g. the auto- 
Backlund transformation, similar to what was done recently for NLS solitons [141 133]. (Scattering 
to zero solutions and criterion for the absence of the MTM solitons were established earlier in 
[181 119] and [34], respectively.) These problems remain open for further studies. 

The paper is organized as follows. Global bounds on the H^ norm of the L^-small solutions are 
obtained in Section 2. In Section 3, we prove that the MTM solitons are non-degenerate minimizers 
of the functional A^ in H^ under the constraints of fixed Q and P for w G {—uq^ojq) with some 
ojo G (0, 1]. Appendix A reports derivation of the conserved quantity R. Appendix B lists a number 
of technical results without proofs. 



2 Global H^ bound on the L^-small solutions 

Thanks to local existence in L^ [8] and the conservation of Q, the L^ solutions are extended globally 
for all t G M with a global bound on the L^ norm of the solution. On the other hand, local existence 
holds also in i7"(M) for any integer n G N but the apriori bounds on the H"^ norm grows at a super- 
exponential rate [161 ^^ . Here we use the conservation of R to obtain the global bound on the H^ 
norm of the solution, which holds for all L^-small solutions. The following theorem gives the main 
result of this section. 

Theorem 2. There is a Qo > such that for all {uq,vo) G H^(W) with ||uo||^2 + ||^o||/,2 ^ Qo, 
there is a positive {uq^vq)- dependent constant C{uq,vo) such that 

\Mh^ + ll^ll^i < C{uo,vo), (8) 

for all t G M. 

Proof. By Sobolev embedding of i/^(M) into L^{M.) for any p >2, the values of Q and R are finite if 
(uo, fo) G H^{M.). To obtain the assertion of the theorem, we need to show that the value of R gives 
an upper bound for the value of ||0^u||^2 + ||c^xw||^2- Then, the standard approximation argument 
in Sobolev space //^(M) yields a conservation of R from the balance equation [see equation (|5U|) in 
Appendix A], whereas the standard continuation argument yields the global bound ([8]). 

The lower bound for R follows from two applications of the Gagliardo-Nirenberg inequality in 
one dimension. For any p > 1, there is a constant Cp > such that 

ll^llii < CpWdMl^'L^^Mllt", u G H\R). (9) 

First, we note that quadratic and sixth-order terms are positive definite. To control the last 
fourth-order terms of R from below, we note that there is a positive constant C (which may change 
from line to another line) such that 



{\u\ + \v\ ){vu + uv)dx 



< C {\\u\\% + \\v\\%) 

< C{\\u\\l, + \\v\\l,){\\d,uh2 + \\d,vh2). 



Because the positive part of it! is quadratic in ||92:^||l2 and ||9a:t;||i2 , the previous bound is sufficient 
to control the last fourth-order terms of R from below. On the other hand, the first fourth-order 
terms like 



{uxU — Uxu)\u\ dx 



< C\\u\\l(i\\dxu\\L2 < C\\u\\l2\\dxu\\l2 



can only be controlled from below if ||w||^2 is sufficiently small. This proves the assertion of the 
theorem. D 

Remark 1. One can try to squeeze the first fourth- order terms of R between the positive quadratic 
and sixth-order terms of R. For example, if reduction v = u is used, these terms of R are estimated 
from below by 



2\\dxu\\ ^2 + 4:\\u\\ j^6 — ^i / {uxU — Uxu)\u\ dx > 2\\dxu\\^2 + 4:\\u\\j^e, — G\\dxu\\i2\ 

Jr 



u 



LB- 



Unfortunately, the lower bound is not positive definite. Therefore, we do not know if the global 
bound |3J can be extended to all (not necessarily L?'-small) solutions of the MTM system |i]j. 



3 H^ orbital stability of solitons 

Critical points of the energy functional H + ujQ with fixed w G (— 1, 1) satisfy the system of first- 
order differential equations 

+i^-uju + v = 2\v\'^u, 



-I 






dx 

The stationary MTM solitons ([2]) correspond to the reduction u = U^ and v = Ui^, where Ui^ is a 
solution of the first-order differential equation 

i-^-LoU + U = 2\U\'^U. (11) 

dx 

Integrating this differential equation with the zero boundary conditions, we obtain MTM solitons 
in the explicit form ([3]). 

Remark 2. Two translational parameters of the MTM solitons ^ are obtained from the gauge 
and space translations: 

u{x,t) ^^ u{x + xo,t)e^°', v{x,t) t-^ v{x + xo,t)e^°', (12) 

where a and xq are real-valued. On the other hand, more general moving MTM solitons must have 
another parameter of velocity c G (— 1, 1), which can be recovered by using the Lorentz transforma- 
tion: 



1/4 



«(x,i)^(l±f) U 



v{x,t)^ U-f 



1/4 



V 



x—ct 
Vl-c2' 


t—cx 
' Vl-c' 


x~ct 


t-cx 
^l-c\ 



(13) 



In what follows, without the loss of generality, we simplify our consideration by working with the 
stationary MTM solitons for c = 0. 

Critical points of the energy functional R + VtQ satisfy the system of second-order differential 
equations 

f ^ + 2i(|n|2 + \v\^)^ + 2iuv^ - 2|t;|2(2|u|2 + \v\^)u + {2\u\^ + \v\^)v + u^v = ^u, 

\ ^-2i{\u\^ + \v\^)^-2iuv^-2\um\u\^ + 2\v\^)v + {\u\^ + 2\v\^)u + v^u = nv. ^ ' 

Using the reduction u = U and v = IJ, we obtain a second-order differential equation 

^ + 6i\U\^^ - 6\U\^U + 3\U\^U + U^ = nU. (15) 

ax^ dx 

Substituting the first-order equation ([TT]) to the second-order equation ([15]) yields the constraint 

(1 - uj'^)u + (2|c/|^ + 2uj\u\'^ -u^ -u^)u = nu, 

which is satisfied by the MTM soliton U = U^^ in the explicit form ([3]) if = 1 — cj^. Therefore, 
the MTM soliton ([3]) is a critical point of the energy functional 

A^ :=i?+(l-w2)Q, wG(-l,l) (16) 



in the energy space i/^(R, C^). 

We shall now prove that there is wq £ (0, 1] such that for any fixed w G (— wq, wq), the critical 
point of A^ is a local non-degenerate minimizer in the constrained space X^^, which is defined as 
an orthogonal complement in i7^(M, C^) of the following complex- valued constraints: 



{u,v) eC^: 
{u,v) eC^: 



{UujU + Uujv) dx = 0, 
[tJ'^u + U'^v) dx = 0, 



(17) 
(18) 



where the prime denotes the derivative of Uui with respect to x. 

The real part of the constraint ()17p is equivalent to the condition that the conserved quantity 
Q is fixed under the perturbation {u,v) to the MTM soliton {Uui,Uuj) at the first order. The 
imaginary part of the constraint (|17p represents the orthogonality of the perturbation (u, v, u, v) to 
the following eigenvector of a linearization operator for the zero eigenvalue, 



9 •" 






(19) 



which is induced by the gauge translation of the MTM soliton {U^^ , U^^ , tJuj , U^^) related to the 
parameter a in the transformation (I12p . 

Similarly, the imaginary part of the constraint (jlSp is equivalent to the condition that the 
conserved quantity P is fixed under the perturbation {u,v) to the MTM soliton {U;^,Uuj) at the 
first order. The real part of the constraint ()18p represents the orthogonality of the perturbation 
{u,v,u,v) to the following eigenvector of a linearization operator for the zero eigenvalue. 



K 



(20) 



which is induced by the space translation of the MTM soliton {Ui^,U^,U^,Ui^) related to the 
parameter xq in the transformation (I12p . 

The following theorem gives the main result of this section. 

Theorem 3. There is coq G (0, 1] such that for any fixed uj G (— a;o,'^o); the Lyapunov functional 
A^ defined by [TS\} is strictly convex at {u,v) = {Uu,,Ui^) in the orthogonal complement of the 
complex-valued constraints pTd and J7gj) in i/^(M, C^). 



To prove Theorem[3l we use a perturbation {u, v, u, v) to the MTM soliton (C4;, f^;, U^, U^) and 
expand the Lyapunov functional A^^ to the quadratic form in {u,v,u,v), which is defined by the 
Hessian operator 

Li 2L2 L2 L3 

2L2 Li L3 L2 

L2 L3 Li 2L2 

L3 L2 2L2 Li 



L 



(21) 



where 



L3 



dx'^ 



4i|C/,. 



dx 



dx 



2f/2 - 2m + 1 



■ OJ 



dUu^ 
dx 



dx 






dx 



U^ - U\ 



Similarly to the case of linearized Dirac equations [9], the 4x4 matrix operator L is diagonalized 
by two 2x2 matrix operators L± by means of the orthogonal similarity transformation 



S'^LS 



L+ 
L_ 







' 1 





-1 







1 





1 





1 


, where S = 














V2 





1 





— 






1 





1 






-1 



The matrix operators L± are found from the reduction of L under the constraints v = ittt: 



L, 



-6a;C/2 






L_ 



2ujW 



(22) 



where 



-^-u\u.\^^ + Q\uX -^ul + wl 



Quj\Uj^ + l 



LJ 



-— -2.|f/.pA 
dx'^ dx 



2 If/,. 



u,t + m-2io\uj' + i 



U! 



and the first-order differential equation (|11|) for Ui^ has been used. Thanks to the exponential decay 
of Uijj to at infinity, by Weyl's Lemma, the continuous spectrum of operators L± is located on the 
semi-infinite interval [1 — a;^,oo) with 1 — o;^ > 0. The following results characterize the discrete 
spectrum of operators L±. 



Proposition 1. For any uj G (—1, 1), we have 

L 



\ ^'" 1 




' ■ 


, L_ 


\ ^^ 1 




" " 


[ul\ 







. -Uu. . 








which represent the eigenvectors M9\) and ^ 
associated with the eigenvectors 



(23) 



In addition, for oj = 0, there is a zero eigenvalue 



cj = : L. 











L_ 



C/o 








(24) 



Proof. Validity of (I23p for any uj G (—1,1) is obtained by direct substitution with the help of 
the differential equations (jlip and (jlSp . Eigenvectors ()19p and ()20p are related to the physical 
symmetries of the MTM system ([1]) with respect to the gauge and space translations. 

Because operators L± are diagonal for a; = 0, the existence of the eigenvectors (|24p follows from 
the existence of the eigenvectors (|23p for w = 0. D 



Lemma 1. For any uj E (—1,1), operator L„ has exactly two eigenvalues below the continuous 
spectrum. Besides the zero eigenvalue associated with the eigenvectors i23\) . L_ also has a positive 
eigenvalue for ui E (0,1) and a negative eigenvalue for uj E (—1,0), which is associated with the 
eigenvector in |^^[ ) for a; = 0. 

Proof. Let us consider the eigenvalue problem L„u = /iu, where u = {u^u) is an eigenvector and 
// is the spectral parameter. Using the transformation 

u{x) = y,(2;)e-*i'o''l^"(^')l''^^' 

where 93 is a new eigenfunction, we obtain an equivalent spectral problem: 



-dl + l-uj"^ - 2w|f/^|2 - 3|C/^|4 
2uj\UJ^ 



-dl + l 






thanks to the fact that 



f^2g2i/;|f/„(x')p 



dx' 



UJ 



I-UJ' 



CO + cosh(2-v/l — uj'^x) 
Because the off-diagonal entries are real, we set 



31^,. 



\U^ 





' V 




" ^ ' 






- /^ 






I ^ J 




[ ^ J 



V'lt := fi^) i ^{x), z := v 1 — Lo^x, fi := {I — a;^)A 

to diagonalize the spectral problem into two uncoupled spectral problems associated with the linear 
Schrodinger operators: 



dz^ 



+ 



3(1-0;^ 



and 



~d^ 



+ 



1 



{oj + cosh(2z))2 

3(1-^2) 



^p+ = X^p+ 



4a; 



(cj + cosh(22;))2 CO + cosh(2z) 
The eigenvector ([23]) in the kernel of L_ yields the eigenfunction 

1 



■0- = Xip- 



(25) 



(26) 



■00 (^) 



(a; + cosh(2z))i/2 



of the spectral problem (|26p for A = 0. Because the eigenfunction -00 is positive definite, the simple 
zero eigenvalue of the spectral problem (j26p is at the bottom of the Schrodinger spectral problem 
for any to E (—1,1), by Sturm's Nodal Theorem (Theorem A in Appendix B). Furthermore, the 

function 

, , . sinh(22;) 

Mz) = —- rTTTT 

UJ + cosh(2z) 

corresponds to the end-point resonance at A = 1 for the spectral problem 



d^ip 



+ 



1 



[1-uj') 



4uj 



{uj + cosh(2z))^ UJ + cosh(2z) 



tp = Xip. 



(27) 



Because the function ^c has exactly one zero, there is only one isolated eigenvalue below the 
continuous spectrum for the spectral problem (j27p (Theorem A in Appendix B). Now the difference 
between the potentials of the spectral problems (pUj) and ([27]) is 



AViz) 



5{l -uj^) 
(a; + cosh(2z))2' 



where AV > for all z G R and uj € (—1,1). By Sturm's Comparison Theorem (Theorem B in 
Appendix B), a solution of the spectral problem (j26p for A = 1, which is bounded as z — )• — oo, has 
exactly one zero. Therefore, the spectral problem (j26p has exactly one isolated eigenvalue for all 
cj G (—1,1) and this is the zero eigenvalue with the eigenfunction ipQ. 

The difference between the potentials of the spectral problems (j25p and (j26p is given by 



AV{z) 



Aoj 



UJ + cosh(2z) 



Since Ay > for u G (0, 1), the spectral problem ()25p has precisely one isolated eigenvalue for 
UJ G (0, 1) (Theorem B in Appendix B) and this eigenvalue is positive (Theorem C in Appendix 
B). On the other hand, since AV < for uo G (—1,0) and ^o > is an eigenfunction of the 
spectral problem (I26p for A = 0, the spectral problem ()25p has at least one negative eigenvalue for 
(^ G (—1,0) (Theorem C in Appendix B). To show that this negative eigenvalue is the only isolated 
eigenvalue of the spectral problem (j25p , we note that 



w + cosh(22;) > a; + 1 + 22;^, z£ 



and consider the spectral problem 

~1^ 



+ 



3(1 -a;2 



(w + l + 2z2)2 



■ip = Xip. 



Rescaling the independent variable z := t-^ y and denoting ip{z) := ip{y), we rewrite 
equivalent form 

dPj: 3 

"d^"(l + y2)2 

It follows that the function 



i-l±il:'|^,(^-')(i + ")^. 



(28) 

in the 

(29) 



V'c(y) 



Vi + y^ 

corresponds to the end-point resonance at A = 1 for the spectral problem 

(fi) 3 



%2 (1 _|_ y 






(30) 



Because the function i/jc has exactly one zero, there is only one isolated eigenvalue below the 
continuous spectrum for the spectral problem (j30p . Because the difference between potentials of 
the spectral problems (j29p and (j30p as well as those of the spectral problems (I25p and (I28p is strictly 
positive for all uj G (—1, 1), by Theorem B in Appendix B, the spectral problem ()25p has exactly 
one isolated eigenvalue for all uj G (—1, 1) and this eigenvalue is negative for uj G (—1,0), zero for 
cj = 0, and positive for cj G (0, 1). D 



Lemma 2. There is coq £ (0,1] such that for any fixed uj G {—ojq^ojq), operator L_|_ has exactly 
two eigenvalues below the continuous spectrum. Besides the zero eigenvalue associated with the 
eigenvector in (EW> ^+ ^^^o has a negative eigenvalue for to G {0,ujq) and a positive eigenvalue for 
UJ G {—ujQ,0), which is associated with the eigenvector in 1^24^ for w = 0. 

Proof. Because the double zero eigenvalue of -L_|_ at a; = is isolated from the continuous spectrum 
located for [l,oo), the assertion of the lemma will follow by the Kato's perturbation theory |22j if 
we can show that the zero eigenvalue is the lowest eigenvalue of L_|_ at w = and the end-point of 
the continuous spectrum does not admit a resonance. 

To develop the perturbation theory, we consider the eigenvalue problem L+u = /xu, where 
u = (n, u) is an eigenvector and /i is the spectral parameter. Using the transformation 

u{x) = ^{x)e-^^!S\^-^^')\'^^' 
where 99 is a new eigenfunction, we obtain an equivalent spectral problem: 



-di + l 



OJ 



- 6w I Uuj 
-6ojW 



SlUu 



-aa + 1 



(jj 



6ojW 

2 6a;|[/.P 



3\U^ 






l^ 



^ 
^ 



where 



W = [/2e6i/o'|C^c.(x')prfx' 

(l +a;cosh(2vT 
= (1--^)^ 



■ Oj'^X 



+ iVl-a;2sinh(2vT 



■ Oj'^X 



(uj + cosh (/ly/l 



LtP'X 



Setting now z := \J\ — u'^x and fi := {1 — uj'^)X, we rewrite the spectral problem in the form 



where 



and 



-92 + i + yi(z) 

V2iz) 



Viiz) :-- 



V2{Z) 



3(1 



CO 



"P 

6uj 



P 
P 



(ci; + cosh(22;))2 (x) + cosh(22;) 



V2iz) :-- 



-6uj 



1 + a; cosh(2z) + i\/l — w^ sinh(22;^ 



{uj + cosh(2z)) 
The eigenvector (j23p in the kernel of L_|_ yields the eigenvector ((^Oi'^o) with 



(31) 



'/5o(^) 



a;sinh(2z) + ivT— ^cosh(22;) 
(tJ + cosh(22;))3/2 



which exists in the spectral problem pip with A = for all to £ (—1, 1)- Now, for w = 0, A = 
is a double zero eigenvalue of the spectral problem (fSTI) . The other eigenvector is {ipq,—(^q) and 
it corresponds to the eigenvector in (|24p. The end-point A = 1 of the continuous spectrum of the 
spectral problem ()3ip does not admit a resonance for oj = 0, which follows from the comparison 
results in Lemma [TJ No other eigenvalues exist for oj = 0. 



10 



To study the splitting of the double zero eigenvalue if cj 7^ 0, we compute the quadratic form 
of the operator on the left-hand side of the spectral problem (j3ip at the vector {ipo, —(^0) to obtain 



2 / {V2 + V2) \ipo\'^dz = -12io 



-3 + 2^2 _^ cosh(4z^ 
(w + cosh(2z))4 



dz. 



Since the integral is positive for a; = 0, Kato's perturbation theory (Theorem D in Appendix B) 
implies that the zero eigenvalue of the spectral problem (j3ip becomes negative for a; < and 
positive for u > with sufficiently small \uj\. D 



Conjecture 1. The spectral problem I131\) has exactly two isolated eigenvalues and no end-point 
resonances for all u) G (—1, 1). The non-zero eigenvalue is positive for all co G (—1,0) and negative 
for all 00 G (0,1). 

Lemma 3. There is ujq G (0, 1] such that for any to G {—ujo,ujq), operators L± have no negative 
eigenvalues and a simple zero eigenvalue in the constrained spaces X± defined by 



X~ 



X. 



u£ L' 



u£ L' 



(UuiU + Uuju) dx = 



{U'^u — U^u) dx = 



(32) 



(33) 



For operator L_, the result extends to all lo G (—1,1). 

Proof. We use Theorem E in Appendix B and compute the value of a in this theorem explicitly 
both for operators L+ and L_. 

For operator L+, the constraint ([32]) yields the vector s = {1/^1, Uu))- By taking derivative of the 
second-order differential equation (jlSp with respect to Vt, we obtain 



hence 



a 



L, 



a,5^ + K*- 



dn 



dn 



dnUo. 



dx 






(34) 



1 d 

2uj doj 



\Ua,\ dx 



2ujy/l 



UJ 



2' 



where we have used the exact expressions $7 = 1 — w^ and ||f/i^|p = arccos(a;). We verify that 
o" > for CO G ( — 1,0) and a < for oo G (0,1). By Lemma O L+ has no negative eigenvalues 
for CO G (— ti;o,0) and has one negative eigenvalue for to G (0,a;o), whereas the eigenvector {U^, 0"^) 
for zero eigenvalue of L+ is orthogonal to the vector s = {Uuj,Uuj). Conditions of Theorem E in 
Appendix B are satisfied and L^ has no negative eigenvalues and a simple zero eigenvalue in the 
constrained space X+ for all co G {—uoq^uoq). Note that the result holds also for a; = 0, since the 
eigenvector {Uq, —Uq) in (p^ does not belong to the constrained space X^ because 



00 



0: / {UoUl, - UoUl,) dx 



(4|C/o| 



r/2 



C/2) dx = -2ij^0. 



(35) 



For operator L_, the constraint (j33p yields the vector s 
equations (llip and ()15p . we obtain 



(C/^, — C/^). By using the differential 



L_ 



1 



-Uu, 



+ 



1 

Aiuo 



U^ 
U^ 



UL 
-u', 



(36) 
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hence 



a = ^\a 



2' "^' Aiuj 



^—{U^Ui-U^U'S\dx 
^ / (4|C/^|4 -Ul-Ul + 4a;|C/^|2) dx 



dx 



l-oj^ I' l+ojcosh{2^/l-uj'^x) 
2oj iiR (w + cosh(2Vl - w22;))2 

2w 

We verify that a < for to G (—1, 0) and cr > for w G (0, 1). By Lemmadl L_ has one negative 
eigenvalue for co G (— a;o,0) and no negative eigenvalues for u G (0,1), whereas the eigenvector 
{Uu,,—Uui) for zero eigenvalue of L_ is orthogonal to the vector s = {U'^,—U'^). Conditions of 
Theorem E in Appendix B are satisfied and L_ has no negative eigenvalues and a simple zero 
eigenvalue in the constrained space X^ for all oj G {—ojq^ojq). Again, the result holds also for 
oj = 0, since the eigenvector {Uq,Uq) of L_ in (|24p does not belong to the constrained space X_ 
because of the same computation (I35p . D 

Remark 3. Solutions of the inhomogeneous equations 1^34^ and i36\) define so-called generalized 
eigenvectors associated with the zero eigenvalue of the spectral stability problem associated with 
MTM solitons of the MTM system ([T]j. These solutions are related to translation of the soliton 
orbit with respect to parameters uj and c. Indeed, from the Lorentz transformation i fi3|] . we realize 
that the solution \3d^) is related to the derivative of the MTM soliton with respect to parameter c at 
c = 0. 

The proof of Theorem [3] follows from Lemma [3] and the fact that the eigenvectors (j23p for the 
zero eigenvalues of L_|_ and L_ are removed by additing additional constraints 



X+ := lue L^(M) : / (C/> + U'^u) dx = L (37) 

X_ := \u(^ L2(]R) : j [U^u - U^u) dx = 1 , (38) 



which are associated with the eigenvectors ()23p . 

Note that the constraints in Xj^ and X^ give real and imaginary parts of the complex-valued 
constraint (J17p , whereas the constraints in X_|_ and X_ give real and imaginary parts of the complex- 
valued constraint (jlSp . Therefore, the Hessian operator L in (12ip is strictly positive under the 
complex- valued constraints (fT7|) and (fTS]) for oj G {—loq^loq) and the proof of Theorem [3] is complete. 

The proof of Theorem[T]is based on the conservation of functionals i?, Q, and P for a solution of 
the massive Thirring model ([T|) in //^(M, C^) and the standard orbital stability arguments (Theorem 
F in Appendix B). 

A Conserved quantities by the inverse scattering method 

The MTM system ([T]) is a compatibility condition of the Lax system 

— <?=L<?, -Qi$=M, (39) 
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where ${x, t) : M x M -;> C^ and L is given bv [231128] : 



iX_ /O 






I,. ,2 I |2x l^ V\ I 



Let us consider a Jost function <j){x;X), which satisfies the boundary condition 



(40) 



hm e-*'=(^)^0(x; A) 



(41) 



where A;(A) := ^(A ^ — A^) € M if A^ € M. This Jost function satisfies the scattering relation as 

X — )• +00, 



(x; A) ~ a(A)e^*^W^ 



+ 6(A)e- 



-ik{X)x 



where a(A) and 6(A) are spectral coefficients for A^ G M. 
Setting 

(j){x; A) 



expiik{X)x+ / x{^''i^)'ix' 



v{x] A) 
with two functions x(x; A) and ^{x] A) satisfying the boundary conditions 

lim x{x]X) = Q and lim z^(x;A) = 



(42) 



(43) 



and taking a limit x — >• oo in (|43p . we obtain 
a(A) = exp ( / x{x]X)dx 



loga(A) 



X{x;\)dx. 



(44) 



The scattering coefficient a(A) does not depend on time f , hence expansion of J_ x{x] X)dx in 
powers of A yields conserved quantities with respect to t [28] . Substituting equation ([^T]) into the 
x-derivative part of the Lax system (|39p , we obtain 



x= 7:(b 



2 I |2 

— \u\ 



)--^(ATJ+ln)., 



where v satisfies a Ricatti equation 



i^^ + i (2A;(A) + 



v\ — \u\ ] u 



>-U 



Au + -n 1 z^^ + 



\/2V A 



r 



Ai; + -li 1 = 



\/2V A 



(45) 



(46) 



To generate two hierarchies of conserved quantities, we consider the formal asymptotic expansion 
of xi?^\ A) in powers and inverse powers of A: 



x(a;;A) = ^A"xn(x), v{x;\) = ^X'vn[x) 



?i=0 



n=l 



and 



xK^ 



oo ^ oo ^ 



(47) 



(48) 



n=0 



n=l 



13 



We set 

/OO /"OO 

Xn{x)dx, I-n ■■= / Xn{x)dx. (49) 

-OO J — OO 

Substitution of the asymptotic expansion of z^ into the Riccati equation (j46p ahows one to determine 
each i^n and z?„ from which equation ([35]) is used to determine Xn and Xn- Let us exphcitly write 
out first conserved quantities 

Io= / {\u\^ + \v\^)dx, 
h = I {—2uxU + ivu + iuv — 2i\u\ \v\ )dx, 
{—2vxV — ivu — iuv + 2i\u\ \v\ )dx, 



/[- 

Jr 



■Aiuuxx — 2{uxV + uvx) + ^u{u\v\ )x + ^Uxu{\u\ +\v\ ) + i{\u\ + \v\ ) 
- 2iuv{\u\'^ + \v\^) - 2ivu{\u\'^ + \v\^) + 4i|'up|'(;|2(|np + \v\^)]dx, 



and 



1-4 = / [Aivvxx — 2{uxV + uvx) + ^v{v\u\ )x + ^Vxv{\u\ +\v\ ) — i{\u\ +\v\ ) 

Jr 

+ 2iuv{\u\'^ + |up) + 2ivu{\u\'^ + \vf) - 4i\u\^ \v\^ {\u\^ + \v\^)]dx. 

We note that /q corresponds to charge Q in ([3]). After integration by parts, I2 + I-2 corresponds 
to momentum P in ([5]) and I2 — I -2 corresponds to Hamiltonian H m. ^. The higher-order 
Hamiltonian R in ([7]) is obtained from I^ — /_4 after integration by parts and dropping the conserved 
quantity Q from the definition of R. 

Using Wolfram's MATHEMATICA, we also obtain the balance equation for R: 

where 

P = \ux\'^ + \vx\'^ - -{uxU-Uxu){\u\'^ + 2\v\'^) + -{vxv -Vxv){2\u\^ + \v\'^) 



and 



-{uv + uv){\u\^ + \v\-^) + 2\u\^\v\^ {\u\ + \v 



\'Ux\^ - \Vx\^ - -{uxU - Uxu){\u\^ + 2\v\^) - -{VxV - Vxv){2\u\^ + \v\^) 
{uv +uv){\u\ —\v\ ). 



14 



B Auxiliary results used in this work 

We are using the following technical results in the main part of this article. In the next four results, 
we consider a linear Schrodinger operator L : H'^{M) — t- L^(M) given by 

L:=-dl + c + V{x), 

where V G L^{R) D L°°(M) and c> is fixed. 

Theorem A |35l Lemma 4.2]. There exists a unique solution of Luq = XqUq for any Xq < c 
such that uq G Hi^^{M.) and limx^^oo e~^'^~'^'^^uo{x) = 1. If uq has n(Ao) zeros on M, then there 
exists exactly n(Ao) eigenvalues of L for any A < Aq. 

Theorem B \35\ Theorem B.IO]. Let u{x; V) be a solution of Lu = cu such that lima;_!._oo u{x; V) 
1. Assume that Vi{x) > V2{x) for all x G M and u{x; V2) has one zero on R. Then, u{x; Vi) has at 
most one zero on M. 

Theorem C |221 Section 1.6.10]. There exists the smallest eigenvalue Aq < c of L if and only if 
Ao = inf / - \(d^uf + a? + Vix)^^'] dx < c. 

u(im{K):\\u\\^2=lJ^2^ J 

In particular, if Aq = for V = Vq and Ay > 0, then Aq ^ for V = Vq± AV . 

Theorem D |22l Section VII. 4. 6]. Let Aq < c be an isolated eigenvalue of L with the eigen- 
function uq G H'^(R). Then, the perturbed operator L := L + IW with AF G L°°(R) has a 
perturbed eigenvalue Aq near Aq and the sign of Aq — Aq coincides with the sign of the quadratic 
form (Ayuo,uo)L2- 

Theorem E |351 Theorem 4.1]. Assume that H \s a, Hilbert space equipped with the inner 
product (•,•). Fix a vector s G -ff. Assume that L is a self-adjoint operator on H such that 
the number of negative eigenvalues of L is n(L), the eigenvectors of L for the zero eigenvalue are 
orthogonal to s, and the rest of the spectrum of L is bounded away from zero. Then, the number 
of negative eigenvalues of L in the constrained space 

Hc:={ueH: (s, u) = 0} 

is defined by the sign of 

a := {L^^s,s). 

If 0" > 0, then the number of negative eigenvalues of L under the constraint is n{L), whereas if 
o" < 0, then this number is n{L) — 1. 

Theorem F [35] Theorem 4.15]. Let X = H^{R,C^) be the energy space for the solution 
Tp := {u,v) of the massive Thirring model ([I]). Let cp^^ := {Uui,Ui^) be a a local non-degenerate 
minimizer of i? in A under the constraints of fixed Q and P for some oj G (—1, 1). Then, 0^^ is 
orbitally stable in X with respect to the time evolution of the MTM system ([T]). In other words, 
for any e > 0, there is (5 > such that if the initial datum satisfies 



then for all t > 0, we have 



inf 11^1,^0 -e*>a.(- + /3)||x<<5. 



inf ||^-e-0^(- + /3)|U<e. 

a,pGK 
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